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Methods for studying droplets in models with quenched disorder are critically examined. Low 
energy excitations in two dimensional models are investigated by finding minimal energy interior 
excitations and by computing the effect of bulk perturbations. The numerical data support the 
assumptions of compact droplets and a single exponent for droplet energy scaling. Analytic calcula- 
tions show how strong corrections to power laws can result when samples and droplets are averaged 
over. Such corrections can explain apparent discrepancies in several previous numerical results for 
spin glasses. 
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Magnets and superconductors are examples of physi- 
cal systems where quenched disorder often plays a dom- 
inant role. Such systems can exhibit hysteresis effects 
and long relaxation times that are the manifestation of 
the large energy barriers created by the quenched disor- 
der. One scenario that makes predictions for the equilib- 
rium and nonequilibrium behavior of disorder dominated 
phases is the droplet or scaling picture |fl]. Predictions 
in this scenario follow from scaling assumptions for the 
energetic and geometric properties of excitations. For 
simple topological reasons, excitations can be defined as 
regions where the configuration is uniformly related by 
a symmetry to a global ground state (e.g., spin flipped 
domains.) In the droplet picture, the low lying excita- 
tions of size I are connected and compact: they have 
volume ~ l df , with dimension df equal to the system di- 
mension d, and the surface to volume ratio decreases as 
I increases. Droplet boundaries are fractal, with surface 
dimension d s < d. The central ansatz is that the prob- 
ability distribution p(A, I) for the energy A of a droplet 
of size I in a given volume ~ l d has a characteristic scale 
~ I 6 . This distribution is argued to have finite weight 
at A = 0. The two exponents 9 and d s can be used, 
for example, to predict many of the properties of a spin 
glass §. This scenario is consistent with numeric results 
for excitations created by modifying boundary conditions 
P,|l. However, other work has suggested that there 
may be more than one important energy scaling expo- 
nent and more complicated geometries for excitations. 
The proposed distinct exponents separately describe (i) 
boundary induced domain wall excitations and (ii) exci- 
tations induced by internal constraints or external fields. 
It has also been suggested that there is distinct scaling 
for large droplets created by modifying the quenched dis- 
order |Q. It is important to understand these claims, as 
they suggest that the standard droplet picture is, at best, 
incomplete. 

To provide perspective, it is useful to investigate in 
detail systems which lend themselves to precise study, 
where some analytic results are known and large systems 



can be simulated efficiently. Results are presented here 
for a 2D elastic medium and a 2D Ising spin glass. Single 
interior droplets, which include a specified central point, 
are computed for the elastic medium. In contrast with 
work on interior droplets in 2D spin glasses Q, a fast, ex- 
act algorithm is used, allowing for precise checks of scal- 
ing. The responses of the elastic medium and the spin 
glass to bulk perturbations are also calculated. The nu- 
merical results for droplet energies and geometrical char- 
acteristics show that logarithmic or small power law cor- 
rections are quite strong. These corrections can be under- 
stood in detail by arguments within the droplet picture. 
Droplets that are not induced by boundary conditions 
are only bounded above by the system size L and below 
by a discretization scale, so that all scales between must 
be considered when computing averages. Corrections to 
scaling for droplets of fixed scale I, such as l^ 1 or L _1 
corrections (e.g., from lattice discreteness) or unknown 
irrelevant operators, might also be considered. However, 
the scale averaging corrections are apparently dominant 
for some quantities. Such corrections lead to an effec- 
tive energy exponent distinct from 9, as boundary condi- 
tion induced domain walls do not have such corrections. 
To remove scale averaging corrections, one can group the 
droplets by scale I and study the geometry and energy as 
a function of / (or I /L if one is interested in large droplets 
||), as L — > oo. With this analysis, the numerical results 
provide strong evidence that the droplets are "compact" , 
with fractal domain walls, and that there is a single en- 
ergy exponent 9. 

One model that I study here is for a two dimensional 
clastic medium, with scalar displacement field u(x), in- 
teracting with quenched periodic disorder. The contin- 
uum energy functional is TL[u(x)\ = J d?x [Vu(x)] 2 + 
V(u(x),x), where V has short range correlations in its 
second argument and is periodic in its first argument, 
V(u(x) + l,x) — V(u(x),x). This model has been used 
for vortex lattices in superconductors, incommensurate 
charge density waves, and crystal growth on a disordered 
substrate [Hi. The continuum model can be discretized 
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on a scale a, where the disorder and elastic energies bal- 
ance H . As an effective degree of freeedom i is pinned 
to a preferred configuration (up to periodic shifts), the 
displacements Ui are of the form rii + /3,-, for integer m 
and fixed Elastic interactions tend to minimize 

nearest neighbor differences in Ui, with excitations of the 
medium being regions displaced relative to the ground 
state. Since the itj are discretized, domain walls separate 
regions relatively shifted by unit amounts. Numerical 
work for zero temperature (T = 0) has determined prop- 
erties of the ground state and the scaling of boundary 
induced domain wall energies 0,[nyn|l|. Finite tem- 
perature simulations [[uj, both Monte Carlo and combi- 
natorial, have shown that the T = phase is stable at 
finite T. This model is thus a useful prototype for models 
with finite T transitions, such as the 3D spin glass |l5|| . 

Another model treated here is the 2D Ising spin glass, 
with Hamiltonian Ji = — X^fe') Jij s i s j> with spins Si — 
±1 on a triangular lattice and Gaussian distributed Jij. 
The ground states {s°} for samples in this model were 
found by a combinatorial method for a standard graph 
representation |0,[L6|]. 

For the elastic medium, minimal energy domain walls 
about the center of a sample were studied on a square 
lattice using a polynomial time algorithm jL7j that cal- 
culated the energy A Q and the droplet boundary. One 
method to characterize the compactness of droplets is to 
compare Ro, the radius of the smallest circle that en- 
closes the droplet, with Ri, the radius of the largest cir- 
cle contained by the boundary vertices. Droplets can be 
studied in spin glasses by finding the ground state and 
then recomputing the ground state {sf} with modified 
couplings J i:j -> J zj - eL- d s° l s° j This bulk per- 

turbation can introduce excitations on all scales. Sample 
excitations are depicted in Fig. |j} 




FIG. 1. (a) Diagram of the geometry of a sample droplet 
in the 2D elastic medium. The droplet is the region which can 
be displaced with minimal energy cost and contains the sam- 
ple center (dot). The radii Ro and Ri are defined in the text, 
(b) Droplets induced by a bulk perturbation (e = 16) in a 64 2 
spin glass sample. The filled areas have spins flipped relative 
to the unperturbed ground state (s° 7^ sf). For the results 
shown in Fig. ^, e w 0.35, 1.3 and droplets rarely intersect. 

The droplet energies are of great interest, as these 
are believed to determine the static correlation func- 
tions at finite temperature and the relaxation to equi- 



librium. Consider the problem of finding the minimal 
energy droplet around the origin in a system of size L 
Q. Assume that there is a factor b which gives a sepa- 
ration of length scales: droplets differing in size by b are 
independent. One independent droplet excitation could 
be excited at each scale, so that sums over all scales must 
be performed (for similar sums, see Refs. p9|JT^ J^|.) At 
each length scale s, s = 1, 2, . . . , log 6 (L), the distribu- 
tion for the energy A (6 s ) has characteristic scale b and 
has finite weight b~ se at A = §. The total density 
of states p(A Q ) for A D < L 6 is then a sum over s, giving 
p(A ) ~ L~ e 1 1 - (b'L) 8 1, for 6 ^ 0, where b' > is set by 
b and the lattice and boundary conditions, which affect 
the lower and upper ends of the sum. The expected min- 
imum value for A Q scales as ~ L e 1 1 — (b'L) e | , the sub- 
dominant term reflecting that the minimal energy droplet 
is chosen from all length scales from 1 to L. The effective 
energy exponent is then (for 9 < 0) 

ff _ dln(A- ) _ 9 

d\n(L) l-(b'Lf K ' 

Applying Eqn. (1) to the 2D spin glass, taking b' = 2 
and L = 16, gives an effective exponent 9 cS — —0.45, 
apparently quite different from the domain wall value 
9 = —0.28 Q and consistent with the alternate energy 
exponent proposed in earlier numerical work Q . The ef- 
fective exponent converges to 9 quite slowly with L (and 
is relatively insensitive to 6'), as 9 is near zero. 

One case where 9 — for domain walls created by 
boundary conditions is the 2D elastic medium. Large 
domain walls can be created by external strains. By sta- 
tistical tilt invariance of the disorder pH] , the change in 
the sample averaged energy can be found by computing 
the elastic energy only, as the change in the sample av- 
eraged pinning energy is zero. Displacing one end of a 
sample by Su = 1 to induce one domain wall gives an 
elastic energy density ~ L~ 2 over the volume L 2 , so that 
the total domain wall energy scales as a constant (9 = 0.) 
This result is consistent with previous numerical simula- 
tions of boundary induced domain walls pfj| . However, 
the mean interior droplet energy A D (L) can be fit over 
a decade with —0.15 < 9 < —0.23, to within a few per- 
cent for smaller L. Arguments similar to those for 9^0 
can be applied to explain this. There are ln(L/a)/ln(6) 
independent scales to choose from, each with identical 
droplet energy distributions (9 = 0). The inverse of the 
minimal droplet energy [A D (L)] _1 is therefore linear in 
ln(L). This result can also be derived using elasticity 
theory. The displacement at the origin of a region of size 
a costs an elastic energy that scales as ~ [ln(L/a)] _1 . By 
tilt symmetry, the pinning can be averaged over, so that 
A Q (L) oc [ln(L/a)] _1 for interior droplets constrained to 
contain the origin. The numerical results are quite con- 
sistent with these expectations, as shown by the two pa- 
rameter fit displayed in Fig. (in addition, the computed 
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probability of generating a droplet of size I is consistent 
with a distribution uniform in ln(Z) [p~7| . ) 

Similar corrections are important for the magnetiza- 
tion m(h) of a spin glass in response to an external field 
h @. For the case 9 ^ and h< 0{L e - d / 2 ), m is found 
by summing over scales the product of the probability 
y^d/i-e f g enera ting a droplet, its expected contribution 
l d / 2 L~ d to the magnetization, and the number of droplets 
(L/l) d at scale I. This gives m tx ft, \{b' h L)~ e — l|, with 
b' h > a constant, to be compared with the uncorrected 
singular piece m cx hL~ e (for numerics, see Ref. ||). The 
size of the corrections are quite similar to those for 9. 
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FIG. 2. Plot of the inverse of the mean droplet energy A 
vs. system size L in the 2D elastic medium, averaged over 
at least 10 4 samples for each L. la error bars are shown. 
The line indicates the fit A D (L) = [0.363 + 0.263 ln(i)] -1 for 
512 > L > 32, with \ 2 = 1-3 for 3 d.o.f. 

The measurement of geometrical quantities, such as 
boundary length and droplet area, can also be strongly 
influenced by scale averaging corrections if one averages 
over all length scales from 1 to L. When the sample 
averaged area A(L) of interior droplets in the 2D elas- 
tic medium is computed as a function of L and the lo- 



cal exponent dj 



d[ln(A(L))]/d[ln(L)} is computed, 
the local dimension is less than two, which might sug- 
gest fractal droplets. This local exponent slowly changes 
with L, though (Fig. The local exponent for A(L) 
is df = 2 - HL/a')]" 1 + 0(L- 2 ), where a' depends 
on the boundary and lattice cutoffs. A useful procedure 
to reduce the corrections to dj S is to separate out the 
scales and plot the droplet area A(Ro,L) as a (binned) 
function of Ro. Changing the order of the averages gives 
local exponents that are much better fit by a constant. 
For L > 32 and R > 8, A{R ) ~ R d J , with bulk droplet 
dimension df — 2.01(2) (Fig. H.) A similar plot for the 
perimeter confirms [|ll] that d s = 1.25(1), with the sur- 
face to volume ratio vanishing as l d ~ d s for large droplets. 
Droplet compactness can also be confirmed by plotting 
the ratio k = Ro / Ri binned according to Ro ; it is 
found that the distribution of k converges, with mean 



(k(R ,L)) = 2.85(5), when L/2 > R > 32. 
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FIG. 3. Estimates of the droplet fractal dimension in a 2D 
elastic medium (la statistical uncertainties shown.) Fixing L, 
binning droplets by linear size Ro , and computing the discrete 
logarithmic derivative to estimate d[ln(j4))]/d[ln(i?o)] gives 
the droplet dimension d^ ft (i?o) values connected by the solid 
lines. From this method df = 2.01(2). Averaging droplet 
area A over samples and computing the discrete logarithmic 
derivative to estimate d[\n(A)]/d[ln(L)] gives the finite size 
estimate d^- (L). The dashed line shows a fit using the cor- 
rection derived in the text with a — 1.7(3). 

Corrections to geometric and energetic quantities are 
also important when computing link overlaps, such as 
those found in comparing the unperturbed Jij ground 
state, {s?}, with the e-perturbed state. The link overlap 
qi is the fraction of link values SiSj on bonds (ij) which 
are unchanged. By summing the contributions over all 
scales (note the small droplets in Fig. |l|(a)), it can be 
shown that the local exponent for the fraction 1 — qi of 
changed bonds behaves as 



d[ln(l- qi )} 
d\n(L) 



= W 



(cLy-^i - 1 • 



(2) 



where /i; = 6 + 2(d — d s ) and c, c' are constants charac- 
terizing the upper and lower cutoffs. The computed local 
exponent for the 2D spin glass is shown in Fig. [| (/x^ ff is 
relatively insensitive to e, at least for 0.35 < e < 1.3.) 
Only for L > 100 does /i° ff approach the large L limit of 
Hi w 1.18 (using the values 9 = —0.28 and d s — 1.27.) 
Similar results are found for the 2D elastic medium [JL7| . 
The exponent f° r spin overlaps q = L~ d ^2 s°s\, with 
1 — q ~ has smaller corrections of this form and 

may well be dominated by corrections to scaling from 
unknown operators or inverse lengths. Numerics |T^ ] 
show that [i converges much more quickly than fii m 
the 2D spin glass. The measurement of Hi (relative to 
\i) has been used by Palassini and Young j7) to conclude 
that a second energy exponent 9' affects the response 
to bulk perturbations in the 3D spin glass. In three di- 
mensions, the scale averaging corrections decrease more 
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quickly with L, as d—fj,i ~ 1.3 compared with the 2D cor- 
rection exponent d — fii~ 0.82, but the system sizes that 
can be simulated are much smaller. It may be that cor- 
rections due to irrelevant variables or possible 1/L effects 
are dominant, but scale averaging corrections clearly con- 
tribute to errors in /i;. A correction of Sfii m —0.2 for 



L ss 8, from similar c and 
elusions of Ref . H . 



would invalidate the con- 
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FIG. 4. Plot of the local exponent /if ff (from the discrete 
derivative between L = 8,16, ... , 256) for the link overlap qi 
in the 2D spin glass as a function of system size. Results 
are shown for distinct e and boundary conditions (open, with 
boundary spins fixed under perturbation, and link periodic.) 
The expected limit at large L is /i; = 1.18(3) (shaded region.) 

In summary, analysis of numerical data provides a pre- 
cise confirmation of the droplet picture in the bulk of a 
sample, both for the scaling of the energies and for the ge- 
ometrical structure of droplets. In comparing the numer- 
ical results with the droplet picture, care must be taken 
to understand where strong corrections might arise. The 
corrections arising from averaging over multiple scales 
can be predicted in some detail and apply to spin glasses 
and other models where 8 is near zero. The corrections 
to averaged quantities such as the droplet energy A D (L) 
are satisfactorily explained, for L 2 > 10 2 systems, by av- 
eraging simple power laws over scales between 1 and L. 
The link overlap, which inherently averages over scales, 
is strongly affected by finite size corrections for small 8 
and d — d s , with effective exponent corrections greater 
than 0.1 for L < 30 in the 2D spin glass. It has been sug- 
gested that, for topological reasons, distinct 9 exponents 
exist only in d > 3 Q. The numerics in this Letter are 
for d = 2, but they and the general analysis suggest that 
large finite size effects strongly affect results in three di- 
mensions. To reduce these types of corrections, data can 
be binned over droplet size I (or over l/L) at fixed system 
size L, checking for convergence by then increasing L. 
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